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How long does it take a random walker to reach a given target point? This quantity, 
known as a first passage time (FPT), has led to a growing number of theoretical inves- 
tigations over the last decade^. The importance of FPTs originates from the crucial 
role played by first encounter properties in various real situations, including transport 
in disordered mediai^i'^, neuron firing dynamics^, spreading of diseases^ or target search 
processes^ii'^'^ . Most methods to determine the FPT properties in confining domains 
have been limited to eflFective ID geometries, or for space dimensions larger than one 
only to homogeneous media^. Here we propose a general theory which allows us to 
accurately evaluate the mean FPT (MFPT) in complex media. Remarkably, this an- 
alytical approach provides a universal scaling dependence of the MFPT on both the 
volume of the confining domain and the source-target distance. This analysis is appli- 
cable to a broad range of stochastic processes characterized by length scale invariant 
properties. Our theoretical predictions are confirmed by numerical simulations for 
several emblematic models of disordered mediai^, fractals^, anomalous diffusion^i and 
scale free networkai^. 

Transport properties are often characterized by the exit time from a sphere toxit, which is the first time a random 
walker reaches any point at a distance r from its starting point. This quantity is well known for Brownian motion 
in euclidean spaces, and has also been evaluated for finitely ramified deterministic fractal o^^i^^ . In these cases, the 
length scale invariant properties of the walker's trajectories play a key role and lead to the scaling form ^cxit oc r''™ , 
which defines the walk dimension dw Interestingly, it has been shown very recently that a large class of complex 
scale free networks are also invariant under a length scale renormalization scheme defined in ref.(— ), even if they are 
of small world type, namely if their diameter scales like the logarithm of the volume. This remarkable property led 
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in particular the authors of ref.(-12.) to characterize the mean exit time in this class of small world networks by a set 
of scaling exponents. 

However, in many situations, the determining quantity is not icxit but rather the FPT of a random walk starting 
from a source point 5 to a given target point T. Indeed the FPT is a key quantity to quantify the kinetics of 
transport limited reaction&iii§ , which encompass not only chemical or biochemical reactionsi^iiS', but also at larger 
scales interactions involving more complex organisms, such as a virus infecting a celU^ or animals searching for food^. 
The relevance of the FPT has also been recently highlighted in ref.(^^) in the context of scale free networks, such 
as social networksSS^ protein interaction networks^^ or metabolic networks^'^. The FPT and the exit time actually 
possess very different properties. Indeed, the exit time is not sensitive to the confinment, since only a sphere of 
radius r is explored by the random walker. On the contrary, an estimation of the time needed to go from one point 
to another, namely the FPT, crucially depends on the confining environment - the MFPT being actually infinite in 
unbounded domains. Here we propose a general theory which provides explicitly the scaling dependence of the MFPT 
on both the volume of the confining domain and the source-target distance. 

Consider a random walker moving in a bounded domain of size N. Let W{r, t\r') be the propagator, i.e. the 
probability density to be at site r at time t, starting from the site r' at time 0, and P(r, t\r') the probability density 
that the first-passage time to reach r, starting from r', is t. These two probability densities are known to be related 
through^^ 



W{rT,t\rs)^ [ P{rT,t'\rs)W{rT,t~t'\rT)dt'. 
Jo 
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After integration over t, this equation gives an exact expression for the MFPT, provided it is finite: 



/rp\ _ H{rT\rT) ~ H{rT\rs) , . 



where 
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H{y\y')^ {W{Y,t\Y')-W,tAY))dt, (3) 

Jo 

and Wstat is the stationary probability distribution (see Supplementary Information (SI) for details). Equation [2] 
is an extension of an analogous form given in ref.(-24), where no quantitative determination of the MFPT could be 
proposed. The main problem at this stage is to determine the unknown function iJ, which is indeed a complicated 
task, since it depends both on the walk's characteristics and on the shape of the domain. A crucial step which allows 
us to go further in the general case is that H turns out to be the pseudo-Green function of the domain^, which in 
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turn is well suited to a quantitative analysis. Indeed, we propose to approximate H by its infinite-space limit, which 
is precisely the usual Green function Gq: 



where Wo is the infinite space propagator (SI). Note that a similar approximation has proven to be satisfactory in 
the standard example of regular diffusion™. We stress that when inserted in equation [21 this form does not lead to a 
severe infinite space approximation of the MFPT, since all the dependence on the domain geometry is now contained 
in the factor 1/Wstat- This approximation is the key step of our derivation and, as we proceed to show, captures 
extremely well the confining effects on MFPTs in complex media. 

We first consider the case of a uniform stationary distribution T4^stat = 1/-^, which is realized as soon as the network 
is undirected and the number of connected neighbors of a node, the degree, is constant. This assumption actually 
underlies many models of transport in complex media, with the notable exception of scale free networks, which will 
be tackled later on. Following ref.(=i), we assume for Wq the standard scaling : 



where the fractal dimension df characterizes the number of sites Nr ~ r ^' within a sphere of radius r and dw has 
been defined previously. This form ensures the normalization of Wq by integration over the whole fractal set. 
A derivation given in SI then yields our central result: 



for r — jr^ — r^l different from 0. Strikingly, the constants A and B do not depend on the confining domain. In 
addition, while A is related to the small scale properties of the walk, we underline that B can be written solely in terms 
of the infinite space scaling function 11 (a precise definition of A and B is given in SI). These expressions therefore 
unveil a universal scaling dependence of the MFPT on the geometrical parameters N and r. 

Several comments are in order. First, we point out that equation © gives the large N asymptotics of the MFPT 
as a function of N and r as independent variables. In particular the volume dependence is linear with N for r fixed 
in any case, which can not be inferred from the standard scaling (T) oc i'*™, L being the characteristic length of 
the domain of order N^/'^i . However, a global rescaling of the problem r ^ Ar, L XL, when applied to equation 
(ini), gives the standard form (T) ^ X^™ for d^ > df and (T) ^ X'^f for d^ < df in accord with ref.(-^'"-). Second, 
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N{A - Br'^^-'-'^f) for d^ < df 
{T) ^ I N{A + Blnr) for d^ = df 
N{A + Br'^-'-'^s) for d^, > df 



(6) 
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equation ([6]) shows two regimes, which rely on infinite space properties of the walk: in the case of compact exploration^ 
{d-uj > df) where each site is eventually visited, the MFPT behaves like (T) cx Nr'^'^'"'^^ ((T) cx iVln?- for c?^ = df) 
at large distance, so that the dependence on the starting point always matters; in the opposite case of non-compact 
exploration (T) tends to a finite value for large r, and the dependence on the starting point is lost. 

We now confirm these analytical results by Monte Carlo simulations and exact enumeration methods applied 
to various models which exemplify the three previous cases, (i) The random barrier model in 2 dimensions'^ is 
a widespread model of transport in disordered systems whose MFPT properties remain widely unexplored. It is 
defined by a lattice random walk with nearest neighbors symmetrical transition rates T distributed according to some 
distribution p(r). Even for a power law distribution p{T) the scaling function n(^) can be shown to be Gaussiarti^ 
{df — d^ ^ 2), which allows us to explicitely compute the constant B and obtain (T) ^ N {A+ {1/27tDcs) Inr). Here 
Dcs is a diffusion constant depending on p{T) which can be determined by an effective medium approximation^*' (SI), 
(ii) The Sierpinski gasket of finite order is a representative example of deterministic fractals, described in Fig. [1] In this 
case^ df = In3/ln2 < In5/ln2 = d^, so that our theory predicts the scaling (T) jVr^'"^"^"^)/'"^ (iu) The Levy 
flight model of anomalous diffusioniii^ is based on a fat-tailed distribution of jump lengths p{l) oc l^^^^ (0 < /? < 2). 
The walk dimension is now d-^j — (3, while the fractal dimension is the dimension of the Euclidian space d . In 
dimensions d > 2, or in ID when /? < 1, one has df > d^, and our theory gives (T) ^ N [A — Br^'"^). 

Figures 2, 3 and 4 reveal an excellent quantitative agreement between the analytical predictions and the numerical 
simulations. Both the volume dependence and the source-target distance dependence are unambiguously captured by 
our theoretical expressions ([S]), as shown by the data collapse of the numerical simulations. We emphasize that the 
very different nature of these examples demonstrates that the range of applicability of our approach, which mainly 
relies on the length scale invariant property of the infinite space propagator ([5]), is wide. 

Remarkably, these analytical results can be extended to scale free networks. The latter are characterized by a 
power-law degree distribution. A wide class of scale free networks has been proven recently to be invariant under a 
length scale renormalization scheme defined in-^: social networks-", the world wide web22., metabolic networks^^, and 
yeast protein interaction networks (PIN)^. Although the standard fractal dimension d/ of these networks is infinite, 
as their diameter scales as IniV, one can show that they are scale invariant in the following sense : they can be covered 
with Nb non overlapping boxes of size Ib with Nb/N ^ Ib'^^ ■ This renormalization property defines an alternative 
scaling exponent called the box dimension c?b, which is actually equal to the fractal dimension defined earlier as long 
as the networks are not of small-world type. A model of scale free ! networks possessing such length scale invariant 
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properties has been defined recursively in ref.(— i^): the network grows by adding m new offspring nodes to each 
existing network node, resulting in well defined modules. In addition, modules are connected to each other through 
X random links (SI). In this case ds — ln(2m + a;)/ln3 and — ln(6m/a; + 3)/ln3. 

For this class of networks, Wstat(r) is not uniform anymore but proportional to the degree fc(r) of the node r. One 
can use the length scale invariant property to infer the following scaling of the infinite space propagator: 



This form, compatible with the symmetry relations proposed in ref(— ), allows us to perform a similar derivation 
which leads for the MFPT to the same result but where df is to be replaced by d^. We applied this formula 
to an example of scale free biological network, the yeast PIN (see figure 1, right), obtained from the filtered yeast 
interactome developed in ref.(— ), and to the model of ref.(— of scale free fractal network. Figure [5] shows that 
this analytical result is in good agreement with numerical simulations on the PIN network. The data collapse over 
various system sizes for the model of scale free fractal networks provides a further validation of our approach, and 
indicates that our theory has a wide range of applications. 
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FIG. 1: Length scale invariant networks. Left: The Sierpinski gasket (here of order 3) is a representa- 
tive example of deterministic fractal. A sample random path from S* to T is shown. Right: The yeast PIN, 
obtained from the filtered yeast interactome developed by (— ). Picture generated by the LaNet-vi software 
(http://xavier.informatics.indiana.edu/lanet-vi/ 1 

FIG. 2: Random barrier model with a transition rate distribution p(r) ~ (a/r)(r/ro)™, with Fq = 1 and a — 0.5. 
The confining domain is a. L x L square with the target point in the middle. Numerical simulations of the MFPT 
rescaled by the volume N, averaged over the disorder, for 3 different domain sizes. The theoretical curve (black line) 
is given by {T)/N ^ {A + (l/27r_Doff) Inr), where the only fitting parameter is A, since DcS is evaluated in SI. 

FIG. 3: Numerical simulations of random walks on a Sierpinski gasket (log/log plot) for 3 different system sizes 
(order 6,7 and 8). For each set of points, the size of the Sierpinski gasket and the target point are fixed (the target 
point correspond to the point T on the Sierpinski gasket of order 3 in fig. [T]) , and the starting point takes various 
positions on the Sierpinski gasket. The plain line corresponds to the theoretical scaling r'^^''^'^^ . 

FIG. 4: Simulations of Levy flights on a 2D square lattice {(3 = 1). The confining domains are a 50 x 50 , 100 x 100 
and 200 x 200 squares, with the target approximately in the middle. The MFPT is presented as a function of the 
source-target distance for different source points. Simulation points are fitted with (T) ^ iV (A — Br^^^). 

FIG. 5: Simulations of random walks on fractal complex networks of small world type. The MFPT on the PIN 
network (plain blue circles) is fitted by (T)/iV ^ (A + Br'^""'^«), with « 2.86 and da « 2.2 as found in^^. We also 
consider 3 examples of the model of networks defined ini^: (m = 3, a; = 1, ds — dw — 1, red symbols and fitting curve), 
(to = 3, a; = 2, (is — d^, = ln(3/2)/ln3, violet symbols and fitting curve) and (m = 3, a; = 3,dB — d^, — 
0, green symbols and fitting curve). For each example, the MFPT (rescaled by the network volume iV = (1 -f x)^ , 
with fc = 3, 4, 5) averaged over the disorder is presented as a function of the source-target distance for different source 
points, and fitted by the theoretical expression {T)/N ^ (yl -|- ij7.i"(3/^)/in3'^_ -y^g gj^j quite surprisingly a scaling 
independent of m. 
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